Sensitivity of the energy levels of singly ionized cobalt to the variation of the fine 
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We use relativistic Hartree-Fock and configuration interaction methods to calculate the depen- 
dence of transition frequencies for singly ionized cobalt on the fine structure constant. The results 
are to be used in the search for variation of the fine structure constant in quasar absorption spectra. 
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I. INTRODUCTION 



Search for variation of fundamental constants is mo- 
tivated by theories unifying gravity with other interac- 
tions as well as by many cosmological models. The search 
spans the whole lifetime of the universe from Big-Bang 
nuclear synthesis to the present-day veryprecise atomic 
clock experiments (see, e.g. reviews No unam- 

biguous manifestation of the variation of fundamental 
constants have been found so far. However, there is large 
amount of data which is consistent with variation of the 
fine structure constant or the ratio of the electron to pro- 
ton mass [l|-|4|). Most of this data comes from the analy- 
sis of the quasar absorption spectra. The analysis of the 
data obtained on the Keck telescope in Hawaii indicate 
that the fine structure constant a might be smaller in 
early universe However, an analysis of the data 

from the VLT telescope in Chile, performed by different 
groups [ll|, [3 gave a null result. There is an intensive 
debate in the literature about possible reasons for the 
disagreement (see. e.g. [H|-[3)- 



The most probable reason for disagreement is the ef- 
fect of some unknown systematics. One of the ways to 
deal with any unknown systematics is to include as many 
atomic lines into analysis as possible and compare the re- 
sults for different lines for consistency. Atomic frequen- 
cies in different atoms and even frequencies of different 
transitions in the same atom depend on the fine struc- 
ture constant very differently. It is extremely unlikely 
that any unknown systematics would behave exactly the 
same way thus mimicking the variation of the fine struc- 
ture constant. 



It was recently brought to our attention that some lines 
of single-ionized cobalt are observed in the quasar absorp- 
tion spectra [3 ■ To include these lines into analysis one 
needs to know how the frequencies of the correspond- 
ing transitions depend on the fine structure constant. 
To reveal this dependence we perform atomic calcula- 
tions following the technique developed in our previous 
works [H [3. 



TABLE I: 
[etc (a.u.)] 



Configurations and effective core polarizabilities 



in the calculations. 


Parity Configuration 


Ctc 


Even 3d 8 


0.4 


Even 3d 7 4s 


0.5177 


Odd 3d 7 4p 


0.448 


Odd 3d 6 4s4p 


0.75 



II. METHOD 



The dependence of atomic frequencies on the fine- 
structure constant a = e 2 /he appears due to relativis- 
tic corrections. In the vicinity of its physical value 
aa = 1/137.036 it is presented in the form 



uj(x) — ujq + qx, 



(1) 



where luq is the present laboratory value of the frequency 
and x = (a/ao) 2 — 1, and q is the coefficient which is to 
be found from atomic calculations. Note that 



dx 



x=0 



To calculate this derivative numerically we use 

oj(x) — uj(—x) 



<1 



2x 



(2) 



(3) 



Here x must be small to exclude non-linear in a 2 terms. 
In the present calculations we use x = 0.01. 

The atomic structure calculations are performed with 
the use of the relativistic Hartree-Fock method (RHF) 
and the configuration interaction technique (CI). The 
RHF self-consistent procedure is done separately for two 
even and two odd configurations presented in Table HI 
The resulting 3d, 4s, and 4p single-electron functions are 
used as a basis for the CI calculations. Note that these 
states are different in different configurations. For ex- 
ample the 3d state in the 3c? 8 configuration is not the 
same as the 3d state in the 3d 7 Ap configuration, etc. See 
Ref. [3 for details. 

The effective Hamiltonian for N v valence electrons has 
the form 



2 
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FIG. 1: Energy levels of cobalt between 45000 cm -1 and 
50000 cm" 1 as functions of (a/ceo) 2 . Dotted line: J = 3, 
short dash line: J = 4, long dash line: J — 5 



FIG. 2: Energy levels of cobalt between 64500 cm -1 and 
67500 cm -1 as functions of (a/ao) 2 . Dotted line: J = 3, 
short dash line: J = 4, long dash line: J — 5 
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(4) 



here, fti(rj) is the one-electron part of the Hamiltonian: 

Ze 2 

hi = ca-p+(P- l)mc 2 h V core + 5V. (5) 

r 

Here a and (3 are Dirac matrixes, V^ore is Hartree-Fock 
potential due to core electrons and 5V is the term which 
simulates the effect of the correlations between core and 
valence electrons. It is often called the polarization po- 
tential and has the form 



SV 



2(r 4 



(0) 



Here a c is the polarization of the core and a is a cut- 
off parameter (we use a = as)- We treat a c as fitting 
parameters and choose their values to reproduce experi- 
mental position of the configurations on the energy scale. 
Corresponding values of a c are presented in Table HI 
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III. RESULTS 

The results of calculations are presented in Table [TT1 
We include experimental and theoretical energy levels, 



FIG. 3: Energy levels of cobalt between 71000 cm and 
75000 cm -1 as functions of (a/ao) 2 . Dotted line: J = 3, 
short dash line: J = 4, long dash line: J — 5 
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TABLE II: Energies (cm" 1 ), p-factors, and (/-coefficients for 
odd states of Co, which are accessible from the ground state 
( 3 F4) via electric dipole transition. 



Experiment" 



Calculations 



State 


Energy 


9 


9nr 


Energy 


9 


<7 


z b n 


45972.17 


1.30 


1.30 


46021 


1.3057 


-872 


2 5 D 3 ° 


47039.27 


1.442 


1.50 


46715 


1.3865 


-444 


z 5 G° 3 


48151.07 


0.927 


0.9167 47221 


0.9498 


50 


z 3 G° 3 


50036.55 


0.80 


0.75 


48481 


0.8572 


368 


z 3 F° 3 


50381.86 


1.055 


1.083 


48148 


1.0910 


-277 


z 3 D° 3 


51512.41* 


1.32 


1.33 


47606 


1.2426 


-1067 


y 5 D§ 


61240.96 


1.505 


1.50 


64844 


1.4971 


-1054 


z 5 P° 3 


63344.50 


1.67 


1.67 


65862 


1.6504 


-977 




63587.01* 


1.35 


1.33 


66075 


1.3407 


-979 


y 3 n 


64360.26 


1.06 


1.08 


65509 


1.0577 


-964 


y 3 G° 3 


65174.89 


0.78 


0.75 


66286 


0.7615 


-280 


z'Ff 


66017.79 


1.01 


1.00 


66678 


1.0111 


-508 


^ 3 D 3 ° 


67524.20* 


1.33 


1.33 


71424 


1.3101 


-1245 


w 3 D° 3 


69060.26 


1.31 


1.333 


73395 


1.3065 


-1026 


x 3 G° 3 


69356.64 


0.78 


0.750 


72125 


0.7578 


-50 


x 3 F° 3 


70457.93 


1.09 


1.083 


73985 


1.0643 


-974 


z 5 F° 4 


45378.85 


1.42 


1.35 


45393 


1.4078 


-1400 


z 5 D° 4 


46320.96 


1.447 


1.50 


45991 


1.4115 


-1100 


z°G^ 


47807.58 


1.154 


1.15 


46732 


1.1564 


-500 


z^G^ 


49348.43* 


1.11 


1.05 


47804 


1.0850 


-300 


z 3 F° 4 


49697.81'' 


1.19 


1.25 


47080 


1.2376 


-900 


y 5 D° 4 


61388.43 


1.51 


1.50 


64773 


1.4953 


-1000 


y 3 n 


63510.40'' 


1.152 


1.25 


65074 


1.2273 


-1311 


z 3 ni 


63792.90* 


0.86 


0.80 


65490 


0.9591 


-957 


Z G4 


64401.46 


1.07 


1.00 


65814 


0.8933 


-522 


y 3 Gl 


65154.18 


1.06 


1.05 


66503 


1.0185 


-290 


x 3 G° 4 


68843.67 


1.05 


1.05 


71639 


1.0439 


-715 


x 3 F° 4 


70186.30* 


1.26 


1.25 


73782 


1.2505 


-962 


y 3 G° 4 


72009.0 


0.70 


0.80 


73090 


0.8322 


-130 


y x G% 


73147.23 


1.04 


1.00 


72734 


0.9816 


-431 


z 5 F° 5 


45197.78 


1.386 


1.40 


45078 


1.3924 


-1568 


z 5 Gl 


47345.94 


1.260 


1.27 


46065 


1.2482 


-1363 


z 3 G° 5 


48556.16'' 


1.194 


1.20 


46872 


1.2248 


-812 


z 3 Ri 


63306.66 


1.03 


1.03 


65002 


1.0312 


-1274 


y 3 G° 5 


64601.72 


1.16 


1.20 


65453 


1.1756 


-1151 




64957.50 


1.05 


1.00 


66126 


1.0212 


-259 


x 3 G% 


68203.39* 


1.18 


1.20 


71242 


1.1985 


-1077 


z 3 ll 


68829.56 


0.81 


0.83 


71767 


0.8396 


-246 



a Ref. [2J 

'States observed in quasar absorption spectra, Ref. 



(/-factors and g-coefficients. Non-relativistic (/-factors are 
given by 

J(J+l)-L(L + l) + 5(5+1) 

9NR = 1 + 2J(jTT) ' (7) 

where J is total momentum of the atom, L is its angular 
momentum and 5 is spin. The g-factors are useful for 
identification of states. As can be seen from Table [TT| ex- 
perimental g-factor are close to the non-relativistic values 
given by ([7]). This justifies non-relativistic notations for 
the states and allows to group them into fine structure 



multiplets. This is why g-factors are often more impor- 
tant then energies for the identification of the states. An 
interesting example of this kind will be discussed below. 

Meanwhile we would like to stress that the accuracy 
for the calculated g-factors can also serve as an indica- 
tor on what kind of accuracy can be expected for the 
g-coefficients. This is because (/-factors are sensitive to 
configuration mixing in a very similar way as the q- 
coefficients [20]. They can even be used to tune config- 
uration mixing to reproduce experimental g-factors and 
therefore improve the accuracy for the g-coefficients pol ] . 
Corresponding procedure can be used when the energies 
of two states with the same total momentum J come very 
close to each other in the vicinity of the physical value of 
a indicating possible level pseudo crossing when energies 
are considered as functions of the ratio (a/ao) 2 - Fig- 
ures 1 1 121 and [3] show energy levels of Co II as functions of 
the fine structure constant from its non-relativistic limit 
a = to the physical value a — ao — 1/137.036. One 
can see that there are many level crossings but most of 
them are at safe distance from the physical point a = ao • 

An interesting example is the mixing of three levels of 
total momentum J = 3: E= 50036.55 cm -1 , £=50381.86 
cm -1 and i?=51512.41 cm -1 . Fig. Q] shows that these 
three levels cross each other at about (a/ao) 2 = 0.5. This 
probably means that mixing should have little effect on 
g-factors or g-coefficients at (a/ao) 2 = 1- However, cal- 
culations do not reproduce correct order of states on the 
energy scale. Here we have to use comparison between 
calculated, experimental and non-relativistic (/-factors to 
decide which fine structure multiplet a calculated state 
should go to. 

In spite of the fact that calculations does not always 
reproduce correct order of states the results for energies 
and g-coefficients are stable with respect of variation of 
the details of the calculation procedure. We estimate 
that the absolute uncertainty in values of g for all states 
is about the same and of the order of 200 cm -1 . There- 
fore, relative uncertainty for smaller g is higher. This 
is because all states belong to the same 3d 7 Ap configu- 
ration and small values of q are due to repulsion from 
lower states of the same total momentum J. Variation 
in the strength of this repulsion or in any other details of 
the calculations are likely to change the values of the g- 
cocfficients for both interacting states by about the same 
amount [2(| • Note however that in cases of poor accuracy 
for g-factors and wrong order of states the uncertainty 
might be slightly higher than 200 cm -1 . 
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